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The dynamic characteristics of annular turbulent seals applied in high-pressure 
turbopumps can be described by stiffness, damping, and inertia coefficients. This 
paper presents an improved procedure for determining these parameters by using 
measurements made with newly developed test equipment. The dynamic system seal, 
consisting of the fluid between the cylindrical surfaces of the rotating shaft and 
the housing, is excited by test forces (input), and the relative motion between the 
surfaces (output) is measured. Transformation of the input and output time signals 
into the frequency domain leads to frequency response functions. An analytical 
model, depending on the seal parameters, is fitted to the measured data in order to 
identify the dynamic coefficients. Some new results are reported that show the 
dependencies of these coefficients with respect to the axial and radial Reynolds 
numbers and the geometrical data of the seal. 


INTRODUCTION 

In high-pressure turbopumps different types of annular seals are used, as 
illustrated in figure 1. The neck or wear ring seals are provided to reduce the 
leakage flow back along the front surface of the impeller face, and the interstage 
seals reduce the leakage from an impeller inlet back along the shaft to the back 
side of the preceding impeller. The geometry of pump seals may be similar to that 
of plain journal bearings, but they have a higher clearance-to-radius ratio and 
normally a fully developed turbulent flow (axial and circumferential) . 

Besides their designed function of reducing the leakage flow, turbulent seals 
have the potential to produce significant forces and have a large influence on the 
bending vibrations of turbopump rotors. To determine more about the pump's dynamic 
behavior, calculation should be carried out in the design stage. For this task the 
machine designer needs to know the dynamic characteristics of seals as expressed by 
inertia, damping, and stiffness coefficients. 

Different theoretical models exist to determine the dynamic seal coefficients 
(refs. 1 to 4), but more experimental data are needed to compare with theoretical 
results' and to modify existing seal models. Until now only a few measured data 
have been available (refs. 5 to 10). The authors have developed a test procedure 
to find the inertia, damping, and stiffness coefficients of annular turbulent seals. 
First measurement results reported in earlier publications (refs. 9, 10) pointed 
out the usefulness of the employed method. Meanwhile a new test rig was built that 
has some improved measurement possibilities and interchangeable seals. However, 
the measurement principle did not change essentially from the first test facility. 


*This research work was supported by "Deutsche Forschungsgemeinschaft”, German 
Federal Republic. 

179 


FtECEDIMQ ?>AQE SLAVIC NOT FiUfgD 



This paper briefly reviews the theoretical seal models and describes the 
employed experimental procedure for the determination of the seal coefficients. 

Some new results are reported, pointing out the dependence of the coefficients on 
the rotational frequency (circumferential Reynolds number), the axial velocity of 
the fluid (axial Reynolds number), and the geometrical data of the seals. They are 
compared with predictions from theoretical models. 


DYNAMIC SEAL COEFFICIENTS - THEORETICAL MODELS 


Most of the existing theoretical seal models have been developed for the sim- 
plest geometry with a fluid surrounded by a cylindrical rotating part (impeller or 
shaft) and a cylindrical housing (fig. 1). Such annular contactless seals separate 
the two spaces with pressure p E and p^, respectively. Caused by the pressure 
difference 4p » pg - p&, the leakage flow in the axial direction is almost always 
turbulent with average velocity V. Because of the rotation of the shaft a 
velocity is superimposed in the circumferential direction. 


From a rotordynamic point of view the various theoretical models usually 
express a relation between the radial forces acting on the rotor and the corre- 
sponding displacements, velocities, and accelerations of the shaft. If the motions 
about a centered position are small, the dynamic system "seal” can be modeled by a 
linear system with stiffness, damping, and inertia coefficients 
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where M is the mass matrix, C the damping matrix, and K the stiffness matrix 
of the seal. 


Black was the first to derive a dynamic model for short annular seals, taking 
into account a leakage relation for flow between concentric rotating cylinders and 
a fully developed shear flow U = RS2/2 (where R is seal radius and Q shaft 
angular velocity) in the circumferential direction but neglecting the pressure- 
induced circumferential flow. His analysis and results are discussed in references 
1 and 2. 

Childs (ref. 3) recently applied Hirs' lubrication equation to the dynamic 
analysis of seals. He used a perturbation method of Hirs' equation to obtain a 
closed-form analytical definition for a short seal. The zeroth-order equations for 
the centered position lead to a leakage pressure drop relation 

Ap = (p E - Pa) ■ (1 + 5 + 2cr) pV 2 / 2 (2) 

with 

u inlet pressure loss coefficient 

a dimensionless friction loss coefficient (L = seal length, (3) 

S s radial seal clearance) , XL/S 
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m (1+m )/2 

X » n R °{1 + 1/Ad] 0 ; b = R /2R (4) 

o a a c 

The friction loss coefficient X depends on the axial Reynolds number R a and the 
circumferential Reynolds number R c , respectively 

R a = 2VS/v; Rc = RQS/v (5) 

defined with the radial clearance S and the fluid viscosity v. The empirical 
coefficients n Q and hiq have to be determined by a test for the special case. 

The first-order perturbation equations lead to a pressure distribution and 
after integration directly to the short seal dynamic coefficients 
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They are in reasonable agreement with prior results of Black (ref. 1) when vq = 0. 
The expressions point out the dependence of the seal coefficients on several param- 
eters, like operating conditions V, S2, the fluid temperature, the friction and 
entry loss coefficients o, X, % with the empirical constants m Q and n Q and 
the seal geometrical data R, L, and S. It can be shown that a few dimensionless 
parameters are sufficient to determine the stiffness, damping, and inertia quanti- 
ties. In the later discussion the two Reynolds numbers R a and R c are used as 
test parameters as well as the geometrical data for seal clearance and seal length. 

In addition to equations (6), the effect of an inlet swirl is also included in 
Child's original paper (ref. 3), pointing out this influence especially to the 
cross-coupled terms. Childs (ref. 4) as well as Black (ref. 2) have also extended 
their analysis to finite-length seals, and Childs has continued his work by deriving 
expressions for seals with different housings and rotor roughnesses. 


IDENTIFICATION OF THE SEAL COEFFICIENTS 

The experimental determination of the inertia, damping, and stiffness coeffi- 
cients of turbulent seals is possible by means of parameter identification. Uorking 
with this method, first of all the dynamic seal model with corresponding mathemati- 
cal equations, particularly frequency response functions, is required. 

In the measurement step test forces (input signals) are applied to the system 
surface in the radial direction, and the relative motions between the two surfaces 
are measured (output signals). From input and output signals mobility frequency 
response functions can be calculated. 

Finally a parameter estimation procedure is applied, requiring a good correla- 
tion between analytical and measured frequency response functions. 


Dynamic Seal Model 

The test rig seal model consists of the fluid, surrounded by the two cylindri- 
cal surfaces of the rotating shaft and the housing, respectively (fig. 2). The 
shaft of the test rig is designed to be very stiff and is supposed to have only one 
degree of freedom concerning the rotation about the axis. On the other hand the 
surface of the housing is movable only in the radial direction, having the two 
degrees of freedom y and z. 

In these circumstances the dynamic system seal can be modeled by means of the 
above presented inertia, damping, and stiffness coefficients. With external time 
dependent forces F y and F z applied to the seal center in the y and z 
directions, the equations of motion are (see eqs. (1) and (6)) 
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( 7 ) 




Corresponding to this two-degree-of-freedom system, a total of four stiffness 
frequency functions exist, assembled in the complex matrix K^in 
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where <■> is the exciter frequency. The inverse functions of equation (8) are 
the mobility frequency response functions 
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represented by the matrix Hy^ n . 


Measurements of Frequency Response Functions 

From the frequency response functions (8) and (9), it is easier to measure the 
mobilities (9). They can be obtained by picking up input and output signals in the 
time domain and transforming them into the frequency domain by fast Fourier trans- 
formation. The ratio of the Fourier- transformed signals leads to the frequency 
response. Figure 3 shows in principle the measurement equipment. The housing is 
excited by a hammer impact, which is equal to an impulse force as broadband excita- 
tion. The relative displacements between the two seal surfaces (housing and shaft) 
are measured with eddy current pickups. After amplification and an analog-to- 
digital conversion the Fourier analyzer calculates the frequency response charac- 
teristics. A bus system transfers the measured data to a digital computer, where 
the unknown seal coefficients are calculated by means of the following identifica- 
tion procedure. 


Estimation of Dynamic Seal Coefficients by Means of Least-Squares 
and Instrumental Variable Method 

For the estimation of the dynamic seal coefficients a least-squares procedure, 
improved by an instrumental variable method, has been applied. From a theoretical 
point of view the product of the mobility matrix H^ n and the stiffness matrix 
Kj c ^ n should be the identity matrix E. An additional error matrix S, caused by 

measurement noise, has to be considered if measured mobilities H$in are combined 
with analytical stiffnessess Kfcin = K - w 2 M + 


SkinScin = I + I 
HfcintK - w 2 M + iwC} = E + S 


( 10 ) 



The complex equation (10) can be rearranged into two real equations 


H^inK - w^H^inM - oiH^inC = E + S r 


H^inK - w^H^inM + wH^inC = 0 + S^ 


( 11 ) 


corresponding to one exciter frequency. In the case of a broadband excitation we 
have as many equations as frequency lines. This results in an overdetermined 
equation system 


A X = E' + S' (12) 

where A consists of the measured frequency response functions related to the 
exciter frequencies <■», X represents the unknown coefficients M, C, and K, E’ 
is a modified identity matrix, and S' is the error matrix. 

Solution of equation (12) in the sense of least squares means to minimize a 
loss function, defined as the Euclidean norm of S* . By deriving the loss function 
with respect to the seal coefficients, the so-called normal equations 

A T A X = A T E' (13) 

are obtained. They can be used to calculate a first set of unknown coefficients 
X 1 , by the matrices M 1 , C 1 , K*, respectively. Figure A shows this branch on the 
left side. However, the use of the described procedure may fail if the noise-to- 
signal ratio in the measurement is too high. 

Fritzen (ref. 11) has suggested an instrumental variable method that is less 
sensitive to noise in the measurement data and reduces the error of the identified 
coefficients. With matrices M 1 , C* , and K 1 from the first step, analytical 
mobility functions can be calculated. Instead of the matrix A T a new matrix W T 
with instrumental variables is built up using the above-mentioned matrices M 1 , C 1 , 
and K 1 . Matrix A T is substituted by W 3 - , and a new set of seal dynamic coeffi- 
cients can be determined with 

W T A X = W T E' (14) 

This procedure is repeated and after each step the actual estimation is compared 
with that of the last step. Measurement and model frequency response functions are 
also compared. The procedure stops if the correlation is satisfactory. Our expe- 
rience is that the combination of the classical least-squares method with the 
instrumental variable method shows much better results than the least-squares 
procedure alone. 


TEST SETUP 

Figure 5 shows a cross section of the actual test rig, which was very similar 
to the one described in references 9 and 10. Two annular seals were integrated 
symmetrically in a very rigid housing. A stiff shaft, driven by an ac motor, 
rotated inside the housing and acted as the second part of the seals. By using 
metal inserts as seal construction parts, we were able to vary the seal length, the 
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radial clearance, and the surface structure. The shaft was fixed to the foundation 
via roller bearings; the housing was connected to this system by two bolt springs 
and eight flexible bars. Other connections between the housing and the shaft were 
two mechanical seals (lip rings) in the end flanges of the casing and the most 
important part, the water film inside the annular seals. 

The hydraulic part of the test setup consisted of a centrifugal pump that fed 
water into the seals. Filter and slide valves regulated the flow and the tempera- 
ture of the water circulation. Flexible hose pipes connected this circuit and the 
movable housing. With this configuration (throughput of the pump, 4.5 m 3 /h; shaft 
speed 6000 rpm) , Reynolds numbers of about 12 000 were reached in the axial direc- 
tion and up to 6500 in the circumferential direction. The water temperature was 
kept constant at 30 °C. 

Besides this mechanical part of the test setup there was the important elec- 
tronic side for the data acquisition and manipulation. Three different groups of 
measured values can be distinguished: 

(1) Data of the fluid 

(2) Data of the exciter 

(3) Response data of the system 

Pressure, temperature, density, and viscosity characterized the fluid state. Sev- 
eral pressure and temperature pickups were distributed over the test plant. The 
fluid velocity determined from the mass flow rate was measured in the supply line, 
and the shaft speed was displayed by the motor control. With these all measured 
data needed to calculate the seal coefficients were available. The excitation 
force, a blow of a hammer, was registered with a piezoelectric accelerometer mounted 
in the hammer head. The third group of data belonged to the motion of the housing 
relative to the shaft. Eddy current pickups surveyed the distance between both 
parts of the seal without contact. Two working planes with two orthogonal measure- 
ment directions permitted the control of housing displacement (translational and 
rotational) . 

The actual data of force and displacement as functions of time were transformed 
to the frequency domain by a digital spectrum analyzer that calculated the transfer 
functions. The parameter identification itself was carried out by desktop computer 
(fig. 3). 


ADAPTATION OF THEORETICAL CONSIDERATIONS TO TEST SYSTEM 

General Remarks 

As shown in equation (7) it is possible to describe the dynamic behavior of 
the seal with differential equations. We can do the same for the test rig, ignoring 
the seal fluid mechanics. With the mass of the housing (m^ = 32.5 kg), the stiff- 
nesses of its attachment (k S y = 10 kN/m; k sz = 14 kN/m) and an assumed damping and 
stiffness matrix for the mechanical seals (k®, c® 1 ), the equation of motion was 
written in the following form: 
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Looking at the complete test apparatus, equation (7) could be replaced with equa- 
tion (15). So we are only able to measure both influences and identify only •’total" 
parameters (fig. 6). Later on we can correct the coefficients (i.e., reduce them 
to the seal coefficients) by using additional information about the system. At the 
moment we deal with the correction terms of the mass and the attachment stiffnesses 
and assume the lip ring coefficients are negligible (some more comments about this 
appear in the last passage of the paper) . 


Realization of Measurements 

Every measurement set consists of four frequency response functions 
because we have a two-degree-of-freedom system. As mentioned 'our input data 
resulted from multiple-impact signals in order to reach a high energy level and a 
good noise-to-signal ratio during the sampling. This force and the system’s 
response were measured simultaneously and later processed in a digital Fourier 
analyzer. The data set obtained in this way represented a single working condition 
of the test system (i.e., fluid velocity, fluid temperature, and shaft speed were 
constant). For all measurements the shaft speed was varied between 2000 and 
6000 rpm, the water volume per time was varied between 2.5 and A. 5 m 3 /h, and the 
temperature was constant at 30 ®C. 


Identification Procedure 

A desktop computer was used for the subsequent treatment of the measured and 
stored data. A program containing both the identification procedure and an in-core 
solver for the determinant equations with respect to the unknown coefficients was 
used to calculate the parameters and the fitted frequency response functions. Both 
functions (measurement and curve fit) could be displayed to estimate the quality of 
the identification algorithm (a typical readout is shown in fig. 7). The next step 
in data processing was to approximate the parameters for a fixed fluid velocity 
with dependence of the shaft speed by using polynomial equations. Only this last 
representation was used in finding the measurement results. 


DISCUSSION OF RESULTS 

The statements written down in this passage refer to figures 8 and 9. In both 
cases the identified values are presented as functions of the Reynolds numbers and 
in relation to the theory in order to compare the different seal geometries. 

Figure 8 demonstrates the influence of growing axial fluid velocity R a and two 
seal clearances (Si = 0.2 mm, S 2 = 0.4 mm); figure 9 demonstrates the parameter 
changes for two seal lengths (Li = 35 mm, L 2 = 23.5 mm) and constant R a . The fol- 
lowing remarks are more general then those figures imply because they are influenced 
by further measurements not given in this paper. 
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Variation of Fluid Velocity and Shaft Speed 


Higher shaft speed led to higher coefficients regardless of geometry. Mass 
and direct damping values were nearly independent of shaft speed, cross-coupled 
stiffness (x stiffness) and cross-coupled damping (x damping) were linearly related 
to shaft speed, and main stiffness was parabolically related. In the case of 
x-stiffness and x-damping the substantial aspect, higher gradients for higher 
velocities, is not as clear as the parallel shifts of the main damping and stiffness 
terms. The general behavior of the values is in good agreement with one of the 
theories (refs. 1 to 4), for both damping values (the x-damping is approximately a 
linear function of the shaft speed, but there is hardly an influence of the axial 
velocity) and the main stiffness coefficient. The inertia term and the x-stiffness 
differ distinctly from the expectations. According to theory the inertia term 
should be only slightly influenced by the velocity, and x-stiffness should show a 
remarkable increase with the flow. 

Comparing the absolute values with the predicted parameters (calculated for 
our test rig with two seals, additional stiffnesses, and the mass of the housing) 
showed that the measurements overestimated the main stiffness (up to 15 percent) 
and the mass (up to 50 percent) and underestimated the x-stiffness clearly (up to 
60 percent) . These errors decreased with larger seal clearances and sometimes 
turned around. Ue may suppose that there could exist an optimal seal geometry with 
"zero error” to the theory. The measured damping values were always in good agree- 
ment with theory. The difference was within a 5 percent range. At this point we 
do not want to evaluate these results, but one remark should be permitted. In the 
case of the mass term there was a disadvantage we had to deal with: its theoreti- 

cal value of about 2 kg for both seals was less then 7 percent of the additional 
mass of the housing. Probably the identification algorithm or the measurement 
system was not sensitive enough to reflect the proper values. 


Influence of Seal Clearance 

The good agreements between measurement and theory are evident for the examples 
in figure 8 for two clearances (S^ =0.2 mm, S£ * 0.4 mm). In all cases the param- 
eters followed the predictions although the results for the x-stiffness and the 
x-damping were better for larger clearances than for smaller one’s. Increasing the 
clearance reduced all coefficients. The main aspects of fluid velocity and shaft 
speed are reproduced for all clearances. 


Influence of Seal Length 

The statements for seal length are the same as in the last paragraph. The 
behavior of all coefficients due to the shaft speed, the fluid velocity, and the 
seal length were comparable to the theory. Sometimes the deviations from the pre- 
calculations were slightly higher, but in general, short seals and larger clearance 
were in better agreement with theory than long seals and small clearances. Probably 
this can be explained by the pressure drop across the seal. Short seals and wide 
clearances produced small pressure drops and so less energy for the system's dyna- 
mic. The result was small stiffnesses and dampings, and so we were able to excite 
the system with the hammer easier as if there was a high pressure drop. Therefore 
the best measurable seal would be a short, wide annular seal. 
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CONCLUSIONS 


From the results we can extract some tendencies, but some results do not fit 
into the theory for many reasons. First, the test system is complex. We only want 
to determine some special characteristics, and the influences of all others should 
be neglected or assumed to be well known in order to find out the right values. 

The lip rings are one great uncertainty. Some additional measurements without seal 
inserts and the exact application of the described method for this special case led 
to stiffness, damping, and inertia terms for the test setup. For example, the mass 
of the housing had been found out with this method with an error less than 5 percent 
to the weighted mass. We found setup damping and stiffness values that showed a 
distinct dependence on the fluid state in the lip rings (i.e., pressure and shaft 
speed clearly influence these parameters). But without searching for the roots of 
these coefficients we corrected the measured seal coefficients with them and found 
some interesting results: 

1. The difference between theory and measurement was smaller than without the 
seal inserts. 

2. The skew symmetry of the matrices improved. 

3. The error in the mass term was constant. 

The lip rings can presently be described as a conservative element with damping 
and stiffness attributes. But this should be verified in the near future. 

The second critical remark touches the inlet swirl. We were not able to meas- 
ure the circumferential velocity of the fluid at the entrance of the seal. For 
short seals we assumed a fully developed Couette flow with an inlet swirl of zero 
because of the geometric configuration of the test rig. This was not the case for 
longer seals, so that could explain once more some differences between theory and 
measurements . 
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Figure 1. - Annular seal types in turbopurrps. 
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Figure 3. - Data acquisition and processing. 
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Figure 4. - Principle of instrumental variable method. 



Figure 5. - Cross section of test rig. 
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Figure 6. - Adaptation of theoretical considerations to real system. 



Figure 7. - Example of measurements and curve-fit data. 
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Figure 8. - Influence of seal clearance S for L = 35 mm and T = 30 °C. 
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Figure 9. - Influence of seal length L for S = 0.2 nun, R = 8424, and T = 30 °C. 
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